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Abstract 

We propose a novel approach to generate chaotic business cycles in a deterministic setting. Rather 
than producing chaos endogenously, we consider aggregate economic models with limit cycles and equilib¬ 
riums, subject them to chaotic exogenous shocks and obtain chaotic cyclical motions. Thus, we emphasize 
that chaotic cycles, which are inevitable in economics, are not only interior properties of economic mod¬ 
els, but also can be considered as a result of interaction of several economical systems. This provides a 
comprehension of chaos (unpredictability, lack of forecasting) and control of chaos as a global economic 
phenomenon from the deterministic point of view. 

We suppose that the results of our paper are contribution to the mixed exogenous-endogenous the¬ 
ories of business cycles in classification by P.A. Samuelson m- Moreover, they demonstrate that the 
irregularity of the extended chaos can be structured, and this distinguishes them from the generalized 
synchronization. The advantage of the knowledge of the structure is that by applying instruments, which 
already have been developed for deterministic chaos one can control the chaos, emphasizing a parameter 
or a type of motion. For the globalization of cyclic chaos phenomenon we utilize new mechanisms such 
that entrainment by chaos, attraction of chaotic cycles by equilibriums and bifurcation of chaotic cycles 
developed in our earlier papers. 

Keywords: Business cycle models; Exogenous shocks; Period-doubling cascade; Attraction of chaotic 
cycles; Chaotic business cycle 

1 Introduction 

Business cycles are a commonly accepted phenomenon in economics. However, we do not actually observe 
perfectly periodic motions in economic variables. Instead, economic data is highly irregular. One way to 
reflect this in economic models is to allow for stochastic processes. Deterministic differential equations 

can also be turned into a better picture of economic reality by introducing chaosQ 

* Corresponding Author Tel.: +90 312 210 5355, Fax: +90 312 210 2972, E-mail: marat@metu.edu.tr 
^There exists a third approach, which is somewhere in-between the two, where iterated function systems generated by 
the optimal policy functions for a class of stochastic growth models converge to invariant distributions with support over 
fractal sets m- 
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Chaotic economic systems can be viewed as unpredictable due to their sensitivity to initial values, 
which makes forecasting extremely difficult na US mi m]. This is known also as the butterfly effect 
m- Devaney |30| proposed that sensitivity in conjunction with other properties, namely transitivity 
and density of periodic solutions, be considered as ingredients of chaos. An alternative way to prove 
the presence of chaos is by observing the period-doubling cascade |38) . This chaos is also sensitive, since 
there are infinitely many solutions with different periods and they are unstable. We utilize these ways of 
observing chaos in our paper. Importantly, irregularity based on theoretical deterministic chaos can be 
visualized in simulations. 

One should remark that it is not only sensitivity that can be considered as a mathematical represen¬ 
tation of unpredictability, but also the existence of infinitely many unstable periodic solutions. Indeed, 
while the presence of a single periodic solution can be accepted as a strong indicator of predictability 
(if one knows the values of the process during the period, then one knows all its future values), with 
infinitely many unstable periodic solutions all the cycles are unstable, and the trajectory of the dynamics 
wanders around, visiting neighborhoods of the cycles in an unpredictable way. That is the reason why in 
the literature the proof of the existence of a periodic-doubling cascade is accepted as evidence of chaos. 
Stabilizing periodic solutions is named in chaos theory as control of chaos. 

Chaos theory could provide a new approach to economic policy-making. Economists believed initially 
that chaotic dynamics is not only unpredictable, but also un-controllable. The results of Ott et ah m 
showed that control of a chaos can be made by a very small corrections of parameters iMiiin]. This and 
related methods have been widely applied to economic models, as exemplified by Holyst et al. m, Kaas 
|49| . Mendes and Mendes [^, Chen and Chen [23] and many others. 

In the classic book m it is observed that while forced oscillator systems naturally emerge in theo¬ 
retical investigations of several technical and physical devices, economic examples for this special family 
of functions have only rarely been provided. The main reason for this deficiency may lie in the fact that 
the necessary periodicity of the dynamic forcing may not be obvious in most economic applications. Our 
proposals are to apply deterministic and chaotic exogenous shocks to economic models and make them 
more realistic. 

One may view chaos (the lack of forecasting) as undesirable in economics, but unavoidable. Hence 
a deterministic economic model is realistic if it exhibits chaotic motions. We suggest considering the 
presence of chaos in a model not only as an indication of its adequacy, but also as a measure of its 
power. Indeed, the presence of chaos implies that the model generates infinitely many aperiodic motions 
and motions with different periods, which are unstable, and consequently easily affected by control and 
sustained in a desirable mode. In other words, deterministic chaos is essential for the flexibility and 
high-speed adjustment of economic models, an indispensable feature in the modern world. 

The principal novelty of our investigation is that we create a chaotic perturbation, plug it in a regular 
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dynamic system, and find that similar chaos is inherited by the solutions of the new system. We call this 
as the input-output mechanism of chaos generation. This approach has been widely applied to differential 
equations before, but for regular inputs. In the studies dull El [HI, the mechanisms for generating 
chaos in systems with asymptotically stable equilibria are provided. In contrast, in [541 1551 [5S1 ED] 
unpredictability in the solutions of differential equations was considered a result of random perturbations 
with small probability. 

P.A. Samuelson 175] accepts purely endogenous theory as “self-generating” cycle. Following this 
opinion we understand chaos as endogenous if it is self-generated by an economic model. One can find 
detailed analysis of the endogenous chaos in books [58l [71188] and paper [T5| , which are very seminal 
sources on the subject. The dynamics arise in duopoly models [78] . in simple ad hoc macroeconomic 
models El ED]. By applying the Li-Yorke theorem it is shown in IHIIT] that an overlapping generations 
model of the Gale type could generate endogenous chaotic cycles. Discrete equations have been applied 
to investigate the presence of chaos in papers ElET], where models representing a capital stock with 
a maximum capital-labor ratio and a Malthusian agrarian economy are investigated. In Ennin] 
endogenous chaotic cycles are demonstrated in growth cycle models. The multiplier-accelerator model 
of Samuelson m has been modified for generation of chaotic endogenous cycles and investigated in 
miMED]- Investigations in Kaldor’s type models, which are originated from mm and finalized in 
m, showed that they could generate endogenous chaos. 

Economists of the first half of the last century already felt a strong need for a theory of irregularities, 
particularly of irregular business cycles. In his classic book, Samuelson [75] observes that while forced 
oscillator systems naturally emerge in theoretical investigations of several technical and physical devices 
and phenomena, economic examples for this special family of functions have only rarely been provided. 
The main reason for this dearth of evidence may lie in the fact that the necessary periodicity of the 
dynamic forcing may not be obvious in most economic applications. That is, economic phenomena do 
not display the kind of regularity that physical phenomena do. Samuelson [75] states that “... in a physical 
system there are grand conservation laws of nature, which guarantee that the system must fall on the thin 
line between stability and instability. But there is nothing in the economic world corresponding to these 
laws...". In a passage Samuelson m suggests that “It is to be stressed that the exogenous impulses which 
keep the cycle alive need not themselves be even quasi-oscillatory in character." Thus, he was already 
talking about irregular business cycles that emerge as a result of irregular exogenous shocks. Moreover, 
he recognized that “most economists are eclectic and prefer a combination of endogenous and exogenous 
theories." Accordingly, in the present paper we consider economic models that admit endogenous business 
cycles and are perturbed by exogenous chaotic disturbances. Examples of models possessing limit cycles 
are Kaldor-Kalecki models and Lienard type equations with relaxation oscillations which are popular in 
economics. Next, the systems are subject to exogenous chaotic disturbances, sensitive and with infinitely 
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many unstable periodic solutions. 

We propose two techniques of obtaining exogenous chaotic cycles as solutions of differential equations. 
In the first approach, an economic model with a limit cycle is perturbed chaotically to produce a chaotic 
business cycle. In the second one, we consider a system with an equilibrium, perturb it by cyclic chaos 
and observe that a chaotic business cycle emerges as a result. While the hrst method is theoretically 
verified in 12], the second method of cyclic chaos generation is new and is demonstrated in our paper 
through simulations. Currently, we study cases where the shocks enter the system additively, but future 
investigations may involve more complex scenarios, where the disturbance enters the main functions of 
the economic model. 

Goodwin |41| argues that the apparent unpredictability of economic systems is due to deterministic 
chaos as much as to exogenous shocks. In this sense, our results can be interpreted as the transmission 
of unpredictability from one economic system to another, and even models that do not admit irregularity 
in isolation can eventually be contaminated with chaos. Thus, we provide support to the idea that 
unpredictability is a global phenomenon in economics, and demonstrate one of the mechanisms for this 
contagion. Considering the current extensive globalisation process, this is a good depiction of reality. 

Our results demonstrate that the control may become not a local (applied to an isolated model) but 
a global phenomenon with strong effectiveness such that control applied to a model, which is realizable 
easily (for example, the logistic map or Feichtinger’s generic model), can be sufficient to rule the process 
in all models joined with the controlled one. Another benefit of our studies is that in literature controls 
are applied to those systems, which are simple and low-dimensional. Control of chaos becomes difficult as 
the dimensions of the systems increase and the construction of Poincare sections becomes complicated. 
Chaos control cannot be achieved if we do not know the period of unstable motion to be controlled. 
In our case, the control is applicable to models of arbitrary dimensions as long as the basic period of 
the generator is known. For these reasons, the possibility to control generated chaos by controlling the 
exogenous shocks that produce the said chaos is appealing. It is especially appealing from a policy¬ 
maker’s point of view, as it offers a cost-effective way to regulate an economic system. 

Control of chaos is nowadays a synonym to the suppression of chaos. Thus our results give another 
way of suppression of chaos. If we find the controllable link (member) in a chain (collection) of connected 
chaotic systems, then we can suppress chaos in the whole chain. This is the effective consequence of our 
studies. 

1.1 Organization of the paper 

The paper is organized in the following way. In the next section we describe the input-output mechanism 
that serves as the basis of chaos extension and formulate two theorems that provide theoretical support 
to the subsequent discussion. In Section |3] economic models with regular motions - stable equilibrium and 
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orbitally stable cycle - are introduced. These models are chaotically perturbed in the following section to 
obtain the main economic dynamics of the paper. More precisely, Section [¥] considers a constellation of 
five economic models connected unilaterally. The extension of chaos near an equilibrium attractor, the 
entrainment by chaos of limit business cycles, the bifurcation of a chaotic cycle, and the attraction of a 
chaotic cycle are the scenarios of the appearance of chaos, and in some cases of chaotic business cycles, 
in economic models that we demonstrate. The effects of applying OGY control Ea to the models will 
also be presented. Section [S] provides simulation results for the entrainment by chaos of limit cycles of 
economic models with time delay. We compare in detail our method of chaos generation with that based 
on the synchronization of chaos |511 iTOl [75] in Section |6| In particular, we argue that chaotic business 
cycles in the paper cannot be obtained by the synchronization of chaotic systems. In Section[7|we discuss 
our results from the point of view of self-organization, and particularly synergetics of Haken [12] ■ We 
summarize the obtained results in the Conclusion. 


2 The Input-Output Mechanism and Its Applications 

To explain the input-output mechanism of chaos generation, let us introduce systems, which we call the 
base-system, the replicator and the generator. They are intensively used in the manuscript. Consider the 
following system of differential equations. 


z' = B{z), (2.1) 

where B : R” —>■ R” is a continuously differentiable function. The system is called the base-system. 

Next, we subdue the base-system to a perturbation, /(<), which will be called an input and obtain 
the following system. 


y' = B{y)+I{t), 


( 2 . 2 ) 


which will be called as the replicator. 

Suppose that the input I{t) admits a certain property, say, it is a bounded function. Assume that 
there exists a unique solution, y{t), of the replicator system (12.2L with the same property. This solution 
is called an output. The process of obtaining the solution y{t) by applying the perturbation I(t) to 
the base-system m is called the input-output mechanism. It is known that for certain base-systems, 
if the input is a periodic, almost periodic, bounded function, then there exists an output that is also 
a periodic, almost periodic, bounded function. In our paper, we consider inputs of a different nature: 
chaotic functions and set of cyclic chaotic functions. The motions that are in the chaotic attractor of the 
Lorenz system m, considered altogether, give us an example of a chaotic set of functions. Each element 
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of this set is considered as a chaotic function. Both a set of functions and a single function can serve as 
an input (as well as an output), and we will use both types of inputs and outputs in this study. 

We consider base-systems of two kinds: (i) systems with asymptotically stable equilibria, (ii) systems 
with limit cycles. In the former case, we will talk about attraction of chaos by equilibria, and in particular, 
attraction of cyclic chaos by equilibria. If the base-system admits a limit cycle, then we talk about the 
entrainment by chaos of limit cycles or just about entrainment by chaos If the limit cycle in a 
base-system is the result of a Hopf bifurcation [33], we will also talk about the bifurcation of the cyclic 
chaos. 

In our previous papers laiaiiiTiisi we analyzed the extension of chaos when the base system possesses 
an asymptotically stable equilibrium. The present paper focuses mostly on the generation of cyclic chaos 
through unilateral coupling of multiple systems. 

The main source of chaos in theory are difference and differential equations. For this reason we 
consider in our manuscript, inputs, which are solutions of some systems of differential or discrete equations 
equations. These systems will be called ^eueratorsy 

Thus, we consider the following system of differential equations, 

x' = G(t, x), (2-3) 

where the function C? : [0, oo) x R"* —>■ is continuous in all of its arguments. We assume that system 

(lO) possesses a chaotic attractor, and we call this system a generator. If x{t) is a solution of the system 
from the chaotic attractor, then we take 


I{t) = eilj{x{t)), 

and use the function I{t) in equation (12.211 . Here, e is a non-zero real number and the function ip : K™ —?> 
R" is continuous. Since we use x{t) as a perturbation in the network (12.2L we call it a chaotic solution. 
The chaotic solutions may be irregular as well as regular (periodic and unstable) [301 ESI [ZZl [ZS] ■ In this 
study we will utilize also the logistic map m as a generator. 

System (ESD is called sensitive if there exist positive numbers cq and A such that for an arbitrary 
positive number 6o and for each chaotic solution x{t) of (12.31) . there exist a chaotic solution x{t) of the 
same system and an interval J C [0, oo), with a length no less than A, such that ||a:(0) — ir(0)|| < Sq and 
||a:(t) — ir(t)|| > eo for all t € J. 

For a given chaotic solution x{t) of (12.31) . let us denote by ?7x(t)(^i yo)i Vo S R", the solution of (12.21) 
with (0, yo) = yo- System (12.21) replicates the sensitivity of (12.31) if there exist positive numbers ei 

^In future work, economic time series that have been tested for the presence of deterministic chaos may be considered 
(see mi [2311^15511551 .) 
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and A such that for an arbitrary positive number (5i and for each solution ?7x(t)(^j2/o)j there exist an 
interval C [0,oo), with a length no less than A, and a solution %(t)(t,yi) such that \\yo — j/i|| < (5i 
and ||?7x(t)(t, 2/o) — %(t)(t,2/i)|| > ei for all t G . Moreover, we say that system (12.21) is chaotic if it 
replicates the sensitivity of (12.31) and the coupled system (12.31) + (12.21) possesses infinitely many unstable 
periodic solutions in a bounded region. 

Next, we will formulate a theorem that forms the mathematical basis of the paper. 

The following conditions are required: 

(Cl) System (j2.1|l admits a non-constant and orbitally stable periodic solution; 

(C2) System (I2.3|l possesses sensitivity and is chaotic through period-doubling cascade; 

(C3) The functions B and G are bounded; 

(C4) There exists a positive number Lb such that 

\\B{zi) - B{z 2 )\\ < Lb ||2;i - Z 2 ||, 

for all zi, Z 2 G K"; 

(C5) There exists a positive number L^ such that 

\\tp{xi) - tp{x 2 )\\ > L^ llxi - X 2 \\ , 

for all xi,X 2 G K™. 

The following assertion is based on the results in [3]. 

Theorem 2.1 If conditions (Cl) —(C5) hold and |e| is sufficiently small, then there exists a neighborhood 
lA of the orbitally stable limit cycle of i2.1\) such that solutions of S2.S\} which start inside lA behave 
chaotically around the limit cycle. That is, the solutions are sensitive and there are infinitely many 
unstable periodic solutions. 

3 Economic Models: The Base Systems 

In what follows, we will require regular systems, that is, models with asymptotically stable equilibria 
or limit cycles, that can be perturbed to generate chaotic business cycles. In this part of the paper we 
propose three economic models to be used as base systems. 
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3.1 Kaldor-Kalecki model with a steady equilibrium 


Consider the following model of an aggregate economy: 


Y' = a[I{Y,K)-S{Y,K)], 

(3.4) 

K' = I{Y, K) - SK, 

where Y is income, K is capital stock, I is gross investment, and S is savings. Income changes pro¬ 
portionally to the excess demand in the goods market, and the second equation is a standard capital 
accumulation equation. The constant depreciation rate 6 and the adjustment coefficient a are positive. 
This model was studied in detail in m and |88] . It admits a stable equilibrium under certain conditions 
on the functions involved. 

Let us consider the following specification of system (limi with /(r, K) = Y-aY^ + bK, S{Y, K) = sY, 


Y' = a[{l - s)r - oY^ bK], 
K' = Y -aY^ + bK- 6K, 


(3.5) 


where the constant parameters satisfy a>0, 6<0, 0<s<l and 0 < 5 < I. 
One can see that a steady state of (USD with positive coordinates 


Y* 


(5(1 — s) + bs 
aS ’ 


s I (5(1 — s) + bs 

'sy ^ 


exists only ii Ss < 6 + bs. 

The transformations Y = y + Y*, K = k + K*, applied to (13.51) . give us the system 


y = a 


, 36s\ ,, „ aS(l — s) + abs , 

2{s-l)- — \y + bk-ay^- 3^ ^- y^ 


,/ /o „ 3bs\ ,, „ aS(l — s) + abs , 

k' = ( 3s - 2- —\y+{b-5)k-ay^-?>\j - - - y^. 


3.2 A model with a business cycle 

We also investigate the idealized macroeconomic model with foreign capital investment, 


S' = aY +pS{k-Y^), 
Y' = v{S + F), 

F' = mS — rY, 


(3.6) 


(3.7) 


where S{t) are savings of households, Y{t) is Gross Domestic Product (GDP), F{t) is foreign capital 
inflow, k is potential GDP, and t is time. If k is set to 1, then Y, S, F are measured as multiples of 
potential output. The parameters represent corresponding ratios: a is the variation of the marginal 
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propensity to save, p is the ratio of capitalised profit, — is the capital-output ratio, m is the capital 

V 

inflow-savings ratio and r is the debt refund-output ratio. The model in (13.711 was introduced by Bouali 
m, and later studied by Bouali et al. m and Pribylova m 
Consider system (lEZl) with specihed coefficients, 

S'' = ay-H0.1S'(l-y2)^ 

F'= 0.5(5'-fF), (3.8) 

F' = 0.195- 0.25F. 

According to [72], the system (13.81) admits Hopf bifurcation at a = oq = 0.25 and an orbitally stable 
cycle appears as a decreases. 

3.3 A Kaldor-Kalecki model 

Let us consider the system. 


Y' = a\l{Y,K)-S{YM. 

K' = I{Y{t-T),K)-5K. (3.9) 

System (13.91) is a Kaldor model with time delay. Kalecki [SD] introduced the idea that there may be 
a time lag between the time an investment decision is made and the time investment is realised. The 
Kaldor-Kalecki model (13.91) was formalised by Krawiec and Szydlowski |55|, where investment depends 
on income at the time investment decisions are taken and on capital stock at the time investment is 
finished. One can find additional information on the models with delay in the papers ISIIHS]. 

We will study the specihcation 


Y' = 1.5[tanh(F) - 0.25A: - (4/3)F], 

K'= tanh(F(t - r)) - O.hK. (3.10) 

According to Zhang and Wei EZ], the model admits an orbitally stable limit cycle for r > 5.4. More 
precisely, the periodic solution appearance follows a Hopf bifurcation so that the origin is asymptotically 
stable if r < 5.4, and the origin loses its stability and the cycle bifurcates from the origin for r > 5.4. 

4 Extension of Chaos in a Constellation of Economical Models 

To provide a comprehensive illustration for the discussion in the previous sections, we will consider a 
constellation of five unilaterally connected economic models denoted by Aj,, fc = 1,..., 5. The topology 
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of the connection is presented in Figure [U and the models are formulated in system (14.111) . We will 
show that the chaos that appears in Ai spreads to all the other models. A 2 serves as a replicator of the 
chaos of Ai and as a generator of chaos in A 3 and A 4 . Model A 4 is a replicator of the chaos of A 2 and a 
generator of chaos in A 3 . 



Figure 1: The connection topology of the systems Ai — A 3 . 


The following is a system of five unidirectionally coupled models Ai — A 3 . 


Kj+i =/rs;j(l-/tj),|^i 

y'l = ( 1 / 8 ) 2 /! - (5/16)/ci - aiyf - 


3 3 -yai 2 


+ vi(t,9), 


K = ( 1 / 4 ) 2 /! - (3/8)/ci - aiyl - 


3v^ 


2 /!, 


> A 2 


2/2 = ( 1 / 3 ) 2/2 - k 2 - 022/2 - + 0-62/!(0 + 


^2 = (l/2)y2 - (5/4)fe2 - 022/1 - 
5" = 0.23T + 0.1S'(l-r2), 

Y' = 0.5{S + F) + 2{yi{t) + 0.5), 

F' = O.IQS'- 0.25y, 

2/3 = (3/5)2/3 - (4/5)fc3 - 032/1 - 


V602 2 
^ 1 — 2 / 2 : 


> A 3 


> A 4 


3 3y^ 2 




yi + O.OlYit), 


k '3 = (7/10)2/3 - (9/10)2=3 - 032/1 - 


3^^ 

vlo 


2/3 > 


> A 3 


(4.11) 


where oi, 02 , 03 are constants and the piecewise constant functions vi(t,9) and 1 = 2 (i,C) defined as 

follows: 


vi{t,9) 


0.019, if 92j <t < 92j+i, 
0 . 002 , if ^ 2 j —! ^ i — ^ 2 j: 


and 


J^2(i,C) 


0.0006, if C 2 j<t<C 2 j+i, 
0.0017, if C 2 j-i<t<C 2 j, 


(4.12) 


(4.13) 


The sequence 9 = {9j} , j € Z, of the discontinuity instants of the function (14.1211 satisfies the relation 
9j = j + Kj, where the sequence {nj} is a solution of the logistic map Ai with kq £ [0,1]. The sequence 
C = {Cj} , j € Z, of the discontinuity instants of (14.131) satisfies the relation Q = 2y/2j for each j. 
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Examples of shocks of the form (14.121) and (14.131) are natural disasters and extreme events in general, 
such as market crashes. They take a finite number of values (an earthquake either happens or not), but 
their timing is irregular or regular. 

4.1 Description of the models Ai to 

Equation Ai is the logistic map, which will be used as the main source of chaos in system (14.111) . The 
interval [0,1] is invariant under the iterations of the map for the parameter values ^ G (0,4], and for 
/r = 3.8 it is chaotic through period-doubling cascade eg. The logistic map plays a very important role in 
many fields of science, particularly in economics. It can be used to describe economic variables. In m the 
logistic map emerges as the law of motion of the price of the non-numeraire good in a simple discrete-time 
model of an exchange economy with two goods under Walrasian tatonnement. Benhabib and Day m 
showed that a logistic map describes optimal consumption in a simple overlapping generations model 
with a quadratic utility function, and Mitra and Sorger m proved that the logistic map can be the 
optimal policy function of a regular dynamic optimisation problem, if and only if the discount factor 
does not exceed 1/16. 

The logistic map Ai is the generator of chaos for the global system (14.111) and as we mentioned 
above, a generator can be not only with continuous dynamics, but also with discrete, and even hybrid, 
i.e., combining both continuous and discrete. In fact the whole model (14.111) is an example of a hybrid 
system. 

System A 2 describes the aggregate economy of Country 1. It is a perturbed Kaldor model (' (13.61) . 
obtained by setting a = 1, s = 1/8, 6 = 1/16 and b = —5/16. In the absence of the perturbation 
function the model possesses an asymptotically stable equilibrium provided that the number 

ui is sufficiently small. One can verify that the associated linear system admits complex conjugate 
eigenvalues (—1 ± i)/ 8 . The function (14.121) describes a rainfall shock that impacts on the agricultural 
sector and through it on the total output. The higher value of vi implies normal rainfall, while the lower 
value is drought, which leads to lower agricultural production and slower output growth. 

Using the results of HEIE], one can state that the chaoticity of the logistic map Ai with /i = 3.8 
makes the function 9) behave chaotically, and system A 2 is chaotic through period-doubling cascade 
for the same value of the parameter That is, it admits infinitely many unstable periodic solutions and 
exhibits sensitivity. For each natural number p, the system possesses an unstable periodic solution with 
period 2p. Next, in its own turn system A 2 is the generator for the systems A^ and A^. 

System A 3 reflects the dynamics of Country 2. It is obtained by using the coefficients a = 1, s = 1/6, 
i5 = 1/4, 5 = — 1 in the Kaldor model (13.61) and by perturbing it with the solutions of A 2 as well as with 
the periodic function (14.131) . The associated linear system has the eigenvalues (—11 ± \/^)/24. In the 
absence of the perturbation terms 0.6yi(t) and 1^2 (t,C) if the number 02 is sufficiently small, the 
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system admits an asymptotically stable equilibrium. The term 0.6yi{t) describes the effect exports from 
Country 2 to Country 1, modelled as a function of the income of Country 1, yi{t), have on the rate of 
change in the income of Country 2. The function V 2 {t, C,) reflects productivity shocks in Country 2, which 
is a binary variable. The higher value of V 2 stands for faster productivity growth, and the lower value 
for slower productivity growth, which leads to slower output growth. 

Since the periodic motions that are embedded in the chaotic attractor of system A 2 with y = 3.8 and 
the function (14.131) have incommensurate periods, one can confirm using the results of [S] that system 
A 3 is chaotic with infinitely many quasi-periodic solutions in the basis. This will be shown through 
simulations in Figure [9] 

System A 4 describes the aggregate economy of Country 3. It is obtained by perturbing system (13.81) 
with the solutions of ^ 2 - It is a replicator with respect to system A 2 , while the term 2{yi{t) + 0.5) is the 
input. This term represents the effect of exports from Country 3 to Country 1, modelled as a function 
of income in Country 1, yl(t), on the rate of growth of income in Country 3. 

In the absence of perturbations, A 4 possesses an orbitally stable limit cycle [71|. Theorem 12. 1 1 implies 
that system A 4 admits chaotic business cycles, provided that the value of the parameter y = 3.8 is used 
in system A 2 . Since the orbitally stable cycle of system (13.81) occurs through a Hopf bifurcation, one can 
talk about the bifurcation of the cyclic chaos. 

System A 3 models the dynamics of Country 4. It is constructed by perturbing the Kaldor model (13.61) 
with the solutions of A 4 , or in economic terms, by perturbing the aggregate economy with exports from 
Country 4 to Country 3, which are a fraction of the income of Country 3, Y{t). The eigenvalues of the 
associated linear system are —1/5 and —I/IO. In the absence of the perturbation term O.OTF(t), the 
system possesses an asymptotically stable equilibrium, for sufficiently small values of 03 . We will make 
use of system A 3 to demonstrate the attraction of chaotic business cycles. 

4.2 Simulations 

In this part of the paper, we will demonstrate numerically the chaotic behavior in system (14.III) . In what 
follows, we will use ai = 3 x I0“®, 02 = I0“®, 03 = 5 x I0“®, y = 3.8 and kq = 0.63. 

Let us start with system A 2 . Setting the initial data yifttf) = 0.12, ki{tf) = 0.08, where to = 0.63, 
we graph in Figure [5] the yi coordinate of system A 2 . It is seen in the figure that system A 2 behaves 
chaotically. 

To show the extension of chaos by system A 3 , we use the solution in Figure [5] as perturbation in system 
A 3 and present in Figure [3] the time series of the 2/2 coordinate of A 3 . The initial data 222 (^ 0 ) = 0.95, 
^ 2 (^ 0 ) = 0.38, where to = 0.63, is used in the simulation. Figure [3] reveals that the chaos of system A 2 
is extended such that the system A 3 also possesses chaos. In order to confirm the extension of chaos 
once more, we depict in Figure H] the projection of the trajectory of the coupled Kaldor system A 2 — A 3 , 
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Figure 2: The graph of the yi coordinate of system A 2 . 


corresponding to the same initial data, on the yi — ki — y 2 space. 



Figure 3: Extension of chaos by system ^ 3 . 



Figure 4: The projection of the chaotic trajectory of the coupled Kaldor-Kalecki system A 2 — A 3 on the 
yi - ki - 2/2 space. 

Next, we continue with system A4. We take into account system A4 with the solution of A2 that is 
represented in Figure^ and show the trajectory of A 4 with S(to) = 1.67, V(to) = 0.94, F(to) = —5.15, 
where tg = 0.63, in Figure 0 One can observe in Figure [S] that the system A 4 admits a chaotic business 
cycle. 

In order to observe the attraction of the cyclic chaos of system A4, we again use the solution of A4 
with S(to) = 1.67, Y{to) = 0.94, F{to) = —5.15, where tg = 0.63, in Ag and depict in Figure | 6 ] the 
trajectory of system Ag with ysitg) = 0.72, kg^tg) = 0.56. It is seen in Figure | 6 ]that the chaotic business 
cycle of A 4 is attracted by Ag, and the cyclic irregular behavior is extended. 
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Figure 5: Chaotic business cycle of system A4. 



Figure 6 : Attraction of cyclic chaos by system A 5 . 


4.3 Control of extended chaos 


The source of the chaotic motions in system (14.111) is the logistic map Ai. Therefore, to control the chaos 
of the entire system, one has to stabilize an unstable periodic solution of the logistic map. The OGY 
control method [ST] is one of the possible ways to do this. We proceed by briefly explaining the method. 

Suppose that the parameter /i in the logistic map Ai is allowed to vary in the range [3.8 — e, 3.8 + e], 
where e is a given small number. That is, it is not possible (say, it is prohibitively costly or practically 
infeasible) to simply shift the value of /r to a level that generates non-chaotic dynamics. Let us consider 
an arbitrary solution {nj} , kq € [ 0 , 1 ], of the map and denote by q = 1,2 ,... ,p, the target unstable 
p—periodic orbit to be stabilized. In the OGY control method eg, at each iteration step j after the 
control mechanism is switched on, we consider the logistic map with the parameter value p = flj, where 


fij = 3.8 




(4.14) 


provided that the number on the right-hand side of the formula (14.1411 belongs to the interval [3.8 — 
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e, 3.8 + e]. In other words, we apply a perturbation in the amount of 


3.8(2k(9) - 1)(kj - 


, . , , to the 

(1 — 

parameter fi = 3.8 of the logistic map, if the trajectory {nj} is sufficiently close to the target peri¬ 
odic orbit. This perturbation makes the map behave regularly so that at each iteration step the orbit 
Kj is forced to be located in a small neighborhood of a previously chosen periodic orbit . Unless 
the parameter perturbation is applied, the orbit kj moves away from due to the instability. If 

> e, we set p,j = 3.8, so that the system evolves at its original parameter 


3.8(2k(J9) - l){Kj - 


k(9)(1 - k(9)) 

value, and wait until the trajectory {nj} enters a sufficiently small neighborhood of the periodic orbit 


<7 = 1 , 2 ,... ,p, such that the inequality —e < 


3.8(2k(«) - 1){kj - 


, ,, , ,, < e holds. If this is the case, 

k(9) (I - «(<?)) 

the control of chaos is not achieved immediately after switching on the control mechanism. Instead, 
there is a transition time before the desired periodic orbit is stabilized. The transition time increases if 
the number e decreases |39| . 

The chaos of system A 2 can be stabilized by controlling an unstable periodic orbit of the logistic map 
Ai, since the map gives rise to the presence of chaos in the system. By applying the OGY control method 
around the fixed point 2.8/3 .8 of the logistic map, we stabilize the corresponding unstable 2—periodic 
solution of system A 2 . The simulation result is seen in Figure [T] We used the same initial data as in 
Figured] It is seen in Figure [7| that the OGY control method successfully controls the chaos of system 
A 2 . The control is switched on at f = ^50 and switched off at t = 0280- The values kq = 0.63 and e = 0.08 
are utilized in the simulation. The control becomes dominant approximately at t = 150 and its effect 
lasts approximately until t = 340, after which the instability becomes dominant and irregular behavior 
develops again. 



Figure 7: The chaos control of system A 2 . The OGY control method is applied around the fixed point 
2.8/3 .8 of the logistic map. The value e = 0.08 is used. 


Next, we will demonstrate the stabilization of an unstable quasi-periodic solution of system A^. We 
suppose that an unstable quasi-periodic solution of A 3 can be stabilized by controlling the chaos of 
system A 2 . We use the solution shown in Figure [7| as the perturbation in system A 3 , and represent in 
Figure [HI the solution of A 3 with y 2 {to) = 0.95, k 2 {to) = 0.38, where to = 0.63. Similarly to system A 2 , 
it seen in the figure that the chaos of A 3 is controlled approximately for 150 < t < 340. 

To reveal that the stabilized solution is indeed quasi-periodic, we depict in Figure [HI the graph of 
the same solution for 200 < t < 300. Figure [HI manifests that application of the OGY control method 
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Figure 8 : The chaos control of system A 3 . It is observable in the figure that controlling the chaos of 
system A 3 makes the chaos of system A 2 to be also controlled. 


to system A 2 makes an unstable quasi-periodic solution of A 3 to be stabilized. On the other hand, the 
stabilized torus of system A 3 is shown in Figure O 



Figure 9: The stabilized quasi-periodic solution of system A3. 


1.148, 



Figure 10: The stabilized torus of system ^ 3 . 


5 Chaotic Business Cycles in Kaldor-Kalecki Model with Time 
Delay 

This section considers the phenomenon of chaos extension by utilizing an economical model with time 
lag (15.1511 . We are devoting a separate discussion to this model, since the result for this case does not 
have theoretical support at the moment. The extension of chaos can be only observed numerically in 
our example, but in the future one could prove the entrainment of the limit cycle by chaos for functional 
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differential equations using the results of the paper [5]. In this section, we will demonstrate numerically 
the formation of chaotic business cycles in the Kaldor-Kalecki model with time delay. 

Let us take into account the system, 


x" + 5{x^ — l)a;' + a; = 5 cos(2.467t),| Bi 

Y' = 1.5 [tanh(F) - 0.25K - (4/3)^] + 0.0045a;(t), (5-15) 

> i?2 

K' = tanh(r(t - r)) - 0.5K. 

Equation Bi is the chaotic Van der Pol oscillator, which is used as the generator system in (15.151) . 
Van der Pol type equations have played a role in economic modelling gnmiiinHi- It is shown by Parlitz 
and Lauterborn [5^ that equation Bi is chaotic through period-doubling cascade. The process of period¬ 
doubling is described by Thompson and Stewart |82| . This implies that there are infinitely many unstable 
periodic solutions of Bi, all with different periods. Due to the absence of stability, any solution that 
starts near the periodic motions behaves irregularly. We will interpret the solution x{t) as an irregular 
productivity shock. 

System B 2 is the Kaldor-Kalecki model and it is the result of the perturbation of the model (13.1011 
of an aggregate economy with a productivity shock. We will observe numerically the appearance of a 
chaotic business cycle, and in particular, the entrainment by chaos of the limit cycle of system (13.1011 . in 
the next simulations. 

Let us take r = 5.5 in i ?2 so that the system possesses an orbitally stable limit cycle in the absence 
of perturbation |57]. We make use of the solution x{t) of Bi with a;(0) = 1.1008, a;'(0) = —1.5546, 
and present in Figure [TT] the solution of B 2 with the initial condition Y{t) = u{t) and K{t) = v{t) for 
t G [—r, 0], where u{t) = —0.057 and v{t) = 0.063 are constant functions. Figure ITT] reveals that the 
dynamics of B 2 exhibits chaotic business cycles. This result shows that our theory of chaotic business 
cycles can be extended to systems with time delay. 
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6 Chaos Extension Versus Synchronization 


Generalized synchronization characterizes the dynamics of a response system that is driven by the output 
of a chaotic driving system [HIMIlHlinilTS]. Suppose that the dynamics of the drive and response are 
governed by the following systems with a skew product structure 

x'=D{x) (6.16) 


and 


y'= R{y,K{x)), (6.17) 

respectively, where x G MP , y G . Synchronization m is said to occur if there exist sets /a,, ly of initial 
conditions and a transformation defined on the chaotic attractor of (16.161) . such that for all x(0) G Ix, 
2 /(0) G ly the relation lim \\y{t) — (j){x{t))\\ = 0 holds. In this case, a motion that starts on Ix x ly 

t—>oo 

collapses onto a manifold M G Ix x ly ot synchronized motions. The transformation (j) is not required 
to exist for the transient trajectories. When (jj is the identity, the identical synchronization takes place 

[2Ql[39]. 

It is formulated by m that generalized synchronization occurs if and only if for all xq G Ix, 2 /i 0 i 2/20 G 
ly, the following asymptotic stability criterion holds: 

lim \\y{t,XQ,yiQ) - y{t,xo,y 2 o)\\ = 0, 

t—¥oo 

where y{t,xo,yio), y{t,Xo,y 2 o) denote the solutions of (I6.17P with y{0,xo,yio) = yio, 2/(0, xq, 2 / 20 ) = 2/20 
and the same x{t), x(0) = xq- 

A numerical method that can be used to investigate coupled systems for generalized synchronization 
is the auxiliary system approach [TJ [32] ■ Let us investigate the coupled economic model A 2 — A 4 for 
generalized synchronization by means of the auxiliary system approach. 

Consider the auxiliary system 


S'q = 0.23Yo + 0.1S'o(l - If), 

Y' = 0.5{So + Fo) + 2y3{t), (6.18) 

= O.lQS'o - 0.25Fo- 

System (16.181) is an identical copy of system A4. 

By marking the trajectory of system A2 — A4 — (16.181) with initial data yi{to) = 0.12, fci(to) = 0.08, 
S{to) = 1.67, Y{to) = 0.94, F{to) = -5.15, S'o(to) = 2.63, Yo{to) = 0.84, Fo(to) = -2.89 at times t = 0j 
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and omitting the first 500 iterations, we obtain the stroboscopic plot whose projection on the Y — Yq 
plane is shown in Figure [T^l Since the plot is not placed on the line Yq = Y, we conclude that generalized 
synchronization does not occur in the couple A 2 — A 4 . 



Figure 12: The auxiliary system approach shows that the systems A 2 and A 4 are not synchronized in 
the generalized sense. 

It is worth noting that generalized synchronization does not take place also in the dynamics of the 
unidirectionally coupled subsystems Bi and B 2 , which are mentioned in Section [SI and this can be 
verified by means of the auxiliary system approach [11 [39] as well. 

7 The Global Unpredictability, Synergetics and Self-Organization 

The idea of the transition of chaos from one system to another, as well as the arrangement of chaos in an 
ordered way, can be viewed through the lens of self-organization |421I65| . Durrenmatt m explained that 
“... a system is self-organizing if it acquires a spatial, temporal or functional structure without specific 
interference from the outside. By ‘specific’ we mean that the structure of functioning is not impressed 
on the system, but the system is acted upon from the outside in a nonspecific fashion." There are three 
approaches to self-organization, namely thermodynamic (dissipative structures), synergetic and the au¬ 
towaves approach. For the theory of dynamical systems (e.g. differential equations) the phenomenon 
means that an autonomous system of equations admits a regular and stable motion (periodic, quasiperi- 
odic, almost periodic). These are what are called autowaves processes [TT] or self-excited oscillations 
[55] in the literature. We are inclined to add to the list one more phenomenon - chaos extension. For 
example, consider the collection of systems Ai, A 2 , ■.., A^ once again, where Ai is the original generator 
of chaos. Because of the connections and the conditions discovered in our analysis, all the other subsys¬ 
tems, A 2 ,. ■., A 5 , are also chaotic. We believe this is a self-organization phenomenon, that is, a coherent 
behavior of a large number of systems [42] . 

Haken [42] . a German theoretical physicist, introduced a new interdisciplinary field of science, syner¬ 
getics, which deals with the origins and the evolution of spatiotemporal structures. Synergetics is based 
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in large part on the dynamical systems theory. One of the crucial features of systems in synergetics is 
self-organization, which was discussed above. According to Haken [32] j the central question in synergetics 
is whether there are general principles which govern the self-organized formation of structures and/or 
functions. The main concepts of the theory are instability, order parameters, and slaving [42] . 

Instability is understood as the formation or collapse of structures (patterns) |5S]. This is very 
common in fluid dynamics, lasers, chemistry and biology |42l 16411651184] . A number of examples of 
instability can be found in the literature on morphogenesis |83j . and pattern formation examples can be 
found in fluid dynamics. The phenomenon is called instability because a former state of fluid transforms 
into a new one, loses its ability to persist, and becomes unstable. One can view the formation of chaos 
in systems A 2 ,..., A 5 in our results as instability. Even though processes in finite dimensional spaces 
are considered, chaotic attractors are assumed to be not single trajectories, but collections of infinitely 
many trajectories with complex topologies. One might say that they are somehow in-between objects 
of ordinary differential equations and partial differential equations. This allows us to also talk about 
dissipative structures m, due to the “density" of the chaotic trajectories in the space. 

Order parameters, when applied to differential equations theory, are those phase variables whose 
behavior produces the main properties of a macroscopic structure and which dominate all other variables 
in the formation, so that the latter can even depend on the order parameters functionally. The dependence 
that is proved (discovered) mathematically is what is called slaving. It is not difficult to see that the 
variables of system Ai are order parameters, and they determine the chaotic behavior of the joined 
systems’ variables. 

8 Conclusion 

We provide examples of models of aggregate economy where the main variables exhibit cycle-like mo¬ 
tion with chaotic elements. Thus, we obtain an irregular business cycle in a deterministic setting. This 
provides a modelling alternative to the business cycle literature relying on stochastic variation in the 
economy. Additionally, our investigation highlights the variety of ways of generating chaos in an eco¬ 
nomic model. Previous work has focused on generating chaos and, in particular, chaotic business cycles 
endogenously (see [111201 [Ml IZIIHE]). Our method of creating chaos has its own relevance for economics, 
since we show the role of exogenous shocks in the appearance of chaos in models that otherwise do not 
exhibit irregular behavior. It can also be said that our work provides a missing link in the research on 
the origins of irregularities in economic time series. While the literature on endogenous chaos was a 
response to the view that exogenous stochastic shocks are the source of fluctuations in the economy (see 
m), this paper is a response to the former, in that it provides a role for exogenous chaotic disturbances 
in producing these fluctuations, and thus completes the circle. 
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Baumol and Benhabib m summarized the significance of chaos research for economics: “Chaos 
theory has at least equal power in providing caveats for both the economic analysis and the policy 
designer. For example, it warns us that apparently random behavior may not be random at all. It 
demonstrates dramatically the dangers of extrapolation and the difficulties that can beset economic 
forecasting generally. It provides the basis for the construction of simple models of the behavior of 
rational agents, showing how even these can yild extremely complex developments. It has served as the 
basis for models of learning behavior and has been shown to arise naturally in a number of standard 
equilibrium models. It offers additional insights about the economic source of oscillations in a number 
of economic models.” 

Indeed, applications of chaos theory have illustrated the possibility of producing complex dynamics 
in deterministic settings Ha [501 ESI mi EHl EH ES], with some papers specifically focusing on building 
“chaotic business cycles " m- Chaos is generated endogenously, and its appearance hinges on the values 
of some crucial parameters of the model. The main novelty of our paper is that we start with a model 
that is not endogenously complex. In one case (model A 4 in the main body), we assume that the system 
has a limit cycle, where the limit cycle is understood to be a closed orbit that is also an attractor m- 
We then subject the model to chaotic exogenous shocks and obtain a perturbed system that admits 
chaotic motions. The chaos emerging around the original limit cycle is cycle-like, and therefore can be 
called a chaotic business cycle. This approach is based on rigorous mathematical theory [HE], and we 
provide numerical simulations. In another case (model in the main body), we subject a system with 
an asymptotically stable equilibrium to chaotic cyclic shocks, which produces a chaotic business cycle in 
the original model, as well. We demonstrate this scenario with simulations, as this approach does not 
have theoretical underpinnings as yet. 

In this paper we show that it is possible to produce a chaotic business cycle in a very natural way 
- take a system of differential equations with a limit cycle as a point of departure, and introduce a 
chaotic exogenous disturbance. An example of an exogenous disturbance is a technology shock to the 
economy which affects output, holding all other variables constant. We describe it using solutions of chaos 
generator models. We use them to demonstrate the proposed approach, and other formulations can be 
studied in future work. For example, one can use actual economic time series, such as commodity prices, 
that have been tested for deterministic chaos [HIIIIEIEH]- Moreover, shocks other than technology 
shocks can be considered, in view of the on-going debate between two literatures supporting and rejecting 
the importance of technology shocks for generating business cycles [HEZlIil]. 

Our results give more theoretical lights on the processes, as we suggest a mathematical apparatus, 
which describe rigorously extension of chaos, increases its complexity, and provides new structures of 
effective control for clusters of economic models. 
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